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Abstract. We show that the symmetry operators for the quantum superintegrable system 
on the 3-sphere with generic 4-parameter potential form a closed quadratic algebra with 6 
linearly independent generators that closes at order 6 (as differential operators). Further 
there is an algebraic relation at order 8 expressing the fact that there are only 5 algebraically 
independent generators. We work out the details of modeling physically relevant irreducible 
representations of the quadratic algebra in terms of divided difference operators in two 
variables. We determine several ON bases for this model including spherical and cylindrical 
bases. These bases are expressed in terms of two variable Wilson and Racah polynomials 
with arbitrary parameters, as defined by Tratnik. The generators for the quadratic algebra 
are expressed in terms of recurrence operators for the one- variable Wilson polynomials. The 
quadratic algebra structure breaks the degeneracy of the space of these polynomials. In an 
earlier paper the authors found a similar characterization of one variable Wilson and Racah 
polynomials in terms of irreducible representations of the quadratic algebra for the quantum 
superintegrable system on the 2-sphere with generic 3-parameter potential. This indicates 
a general relationship between 2nd order superintegrable systems and discrete orthogonal 
polynomials. 

Key words: supcrintcgrability; quadratic algebras; multivariable Wilson polynomials; mul- 
tivariable Racah polynomials 
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1 Introduction 

We define an n-dimensional classical superintegrable system to be an integrable Hamiltonian 
system that not only possesses n mutually Poisson - commuting constants of the motion, but in 
addition, the Hamiltonian Poisson-commutes with 2n — 1 functions on the phase space that 
are globally defined and polynomial in the momenta. Similarly, we define a quantum su- 
perintegrable system to be a quantum Hamiltonian which is one of a set of n algebraically 

*This paper is a contribution to the Special Issue "Relationship of Orthogonal Polynomials and Spe- 
cial Functions with Quantum Groups and Integrable Systems". The full collection is available at 
|http://www.emis.de/journals/SIGMA/OPSF.html| 
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independent mutually commuting differential operators, and that commutes with a set of 
2n — 1 independent differential operators of finite order. We restrict to classical systems of 

n 

the form % = Yl 9^PiPj + ^ an d quantum systems H = A n + V. These systems, inclu- 

ding the classical Kepler [1] and anisotropic oscillator systems and the quantum anisotropic 
oscillator and hydrogen atom have great historical importance, due to their remarkable proper- 
ties [21 El HI [5, 6j. One modern practical application among many is the Hohmann transfer, 
a fundamental tool for the positioning of earth satellites and for celestial navigation in gene- 
ral, which is based on the superintegrability of the Kepler system [7]. The order of a clas- 
sical superintegrable system is the maximum order of the generating constants of the motion 
(with the Hamiltonian excluded) as a polynomial in the momenta, and the order of a quan- 
tum superintegrable system is the maximum order of the quantum symmetries as differential 
operators. 

Systems of 2nd order have been well studied and there is now a structure and classification 
theory [H El [101 CEH ESI [13], especially for the cases n = 2, 3. For 3rd and higher order superin- 
tegrable systems there have been recent dramatic advances but no structure and classification 
theory asyet[ia[l5l[l6l[I71[l8l[ig[^ 

The potential V corresponding to a 2nd order superintegrable system, classical or quan- 
tum, on an n-dimensional conformally flat manifold depends linearly on several parameters 
in general and can be shown to generate a vector space of dimension < n + 2. (One di- 
mension corresponds to the trivial addition of a constant to the potential and usually isn't 
included in a parameter count.) If the maximum is achieved, the potential is called non- 
degenerate. There is an invertible mapping between superintegrable systems on different 
manifolds, called the Stackel transform, which preserves the structure of the algebra gene- 
rated by the symmetries. In the cases n = 2,3 it is known that all nondegenerate 2nd 
order superintegrable systems are Stackel equivalent to a system on a constant curvature 
space |30[ 13 lj . An important fact for 2D systems is that all systems can be obtained from 
one generic superintegrable system on the complex 2-sphere by appropriately chosen limit pro- 
cesses, e.g. \32\ 133] . The use of these processes in separation of variables methods for wave 
and Helmholtz equations in n dimensions was pioneered by Bocher [34] . For n = 3 it ap- 
pears that all nondegenerate 3D systems can be obtained from one generic superintegrable 
system on the complex 3-sphere by similar limiting processes, but the proof is not yet com- 
plete [hies]. 

For n = 2 we define the generic sphere system by embedding of the unit 2-sphere xl+x^+x 2 = 
1 in three dimensional flat space. Then the Hamiltonian operator is 

3 

H = ^2 (xidj - Xjdif + ^2 — > di = d Xi . 

1<«<7<3 k=l °° k 

The 3 operators that generate the symmetries are L\ = Li%, L2 = £13, £3 = -^23 where 
Lij = Lji — [x^dj XjOj^ -\- 2 ~t~ 



X? X 2 a 



for 1 < i < j < 4. Here 

3 



H = E L ij + ^a k = H + V, V=4 + 4 + 4- 

1<«<J<3 k=l 12 3 



From the general structure theory for 2D 2nd order superintegrable systems with nonde- 
generate potential we know that the 3 defining symmetries will generate a symmetry algebra 
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(a quadratic algebra) by taking operator commutators, which closes at order 6, [36]. That is, all 
possible symmetries can be written as symmetrized operator polynomials in the basis generators 
and in the 3rd order commutator R, where R occurs at most linearly. In particular, the dimen- 
sion of the space of truly 2nd order symmetries for the Hamiltonian operator is 3, for the 3rd 
order symmetries it is 1, for the 4th order symmetries it is 6, and for the 6th order symmetries 
it is 10. For the generic 2-sphere quantum system the structure equations can be put in the 
symmetric form [12] 

e ijk [Li,R} = 4:{Li,L k } - 4{Li,Lj} - (8 + 16^)^ + (8 + 16a k )L k + 8(oj - a k ), (1.1) 

R 2 = -{L U L2, L 3 } - (16ai + 12)L\ - (16a 2 + 12)L 2 2 - (16a 3 + 12)L\ 

52 11 
+ — ({L 1: L 2 } + {L 2 , L 3 } + {L 3 , la}) + -(16 + 176^)^ + -(16 + 176a 2 )L 2 

1 32 
+ -(16 + 176a 3 )L 3 + — (ai + a 2 + a 3 ) + 48(aia 2 + a 2 a 3 + a 3 ai) + 64aia 2 a 3 . (1.2) 



Here is the pure skew-symmetric tensor, R = [Li,L 2 ] and {Li,Lj} = LiLj + LjLi with 
an analogous definition of {Li,L 2 ,L 3 } as a symmetrized sum of 6 terms. In practice we will 
substitute L 3 = H — L\ — L 2 — a\ — a 2 — a 3 into these equations. 

In [12] we started from first principles and worked out some families of finite and 
infinite dimensional irreducible representations of the quadratic algebra with structure rela- 
tions f 1 1 . 1 p . (jl.2p . including those that corresponded to the bounded states of the associated 
quantum mechanical problem on the 2-sphere. Then we found 1-variable models of these repre- 
sentations in which the generators Li acted as divided difference operators in the variable t on 
a space of polynomials in t 2 . The eigenfunctions of one of the operators Lj turned out to be the 
Wilson and Racah polynomials in their full generality. In essence, this described an isomorphism 
between the quadratic algebra of the generic quantum superintegrable system on the 2-sphere 
and the quadratic algebra generated by the Wilson polynomials. 

The present paper is concerned with the extension of these results to the 3-sphere, where 
the situation is much more complicated. From the general structure theory for 3D 2nd or- 
der superintegrable systems with nondegenerate potential we know that although there are 
2n — 1 = 5 algebraically independent 2nd order generators, there must exist a 6th 2nd or- 
der symmetry such that the 6 symmetries are linearly independent and generate a quadratic 
algebra that closes at order 6 [37] , (We call this the 5 ==> 6 Theorem.) Thus, all pos- 
sible symmetries can be written as symmetrized operator polynomials in the basis genera- 
tors and in the four 3rd order commutators where the R{ occur at most linearly. In 
particular, the dimension of the space of truly 2nd order symmetries is 3, for the 3rd or- 
der symmetries is 4, for the 4th order symmetries it is 21, and for the 6th order sym- 
metries it is 56. In 3D there are 5 algebraically independent, but 6 linearly independent, 
generators. The algebra again closes at 6th order, but in addition there is an identity at 
8th order that relates the 6 algebraically dependent generators. The representation the- 
ory of such quadratic algebras is much more complicated and we work out a very impor- 
tant instance of it here. In this case we will find an intimate relationship between these 
representations and Tratnik's 2-variable Wilson and Racah polynomials in their full genera- 
lity PHIEEIHD]. 

For nD nondegenerate systems there are 2n — 1 functionally independent but n(n + l)/2 
linearly independent generators for the quadratic algebra. We expect that the relationships 
developed here will extend to n-spheres although the results will be of increasing comple- 
xity. 
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2 The quantum superintegrable system on the 3-sphere 

We define the Hamiltonian operator via the embedding of the unit 3-sphere xf + x^ + x^ + x^ = 1 
in four-dimensional flat space 

4 

H= (^d j -x j d i ) 2 + ^2^, di = d Xi . (2.1) 

l<i<j<4 k=l Xk 

A basis for the second order constants of the motion is 

CLX^ ? 
L^j — Lji — {x^dj Xjdi) -\- 2 ~t~ 2 ' 

X i Xj 

for 1 < i < j < 4. Here 

4 

l<i<j<4 k=l 

In the following i, j, k, I are pairwise distinct integers such that 1 < i,j,k,£ < 4, and 
is the completely skew-symmetric tensor such that = 1 if i < j < k. There are 4 linearly 
independent commutators of the second order symmetries (no sum on repeated indices): 

Re = £ijk[Lij, Ljk]. 

This implies, for example, that 

Ri = [-^23^34] = —[-^24,-^34] = —[^235-^24]- 

Also 

[Lij,L k i] = 0. 

Here we define the commutator of linear operators F, G by [F, G] = FG — GF. The structure 
equations can be worked out via a relatively straightforward but tedious process. We get the 
following results. 

The fourth order structure equations are 

[Lij,Rj] = 4ei£k{{Lik, Lji} — {La,Lj k } + La — Lit + Ljk — Lji), 

[Lij,Rk] = 4:eije({Lij,L ie - L je } + (2 + Aaj)L ie - (2 + 4ai)L je + 2a { - 2aj). 

Here {F, G} = FG + GF. 

The fifth order structure equations (obtainable directly from the fourth order equations and 
the Jacobi identity) are 

[Re, R k ] = 4,e iki (Ri - {L lj: Ri}) + Ae jM (Rj - {Lij,Rj}). 

The sixth order structure equations are 

Rj = ^{Lij , L ik , L jk } - (12 + 16a k )Ll - (12 + 16 ai )L% - (12 + 16^)^ 

52 /16 176 \ /16 176 \ 
+ -^-{{Lij, L ik + Lj k \ + {L ik , Lj k \) + I — + —^~ a k I Lij + I — -| —di I Lj k 

/16 176 \ _ io/ , 32. 

+ I — + "^""i J L ij + 64ajaja fc + 48(ajOj + Oj-a fc + afcaj) + — (a^ + aj + a fc ), 
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lki ^ kl {R h Rj} — -({Lu,Lj k , L ke } + {L ik , Lj£, L ki } — {Lij,L ki , L ki }) + — {L ik , Lji} 

26 44 2 

+ -^-{Lie, Lj k } + —{Lij, L k i\ + 4L M — 2{Ljg + Lj k + La + L ik , L k g\ 

32 

- (6 + 8ag){L ik , L jk } - (6 + 8a k ){L ie , L je } - — L M 



+ ( — + 24a fc + 24a^ + 32a k a e ) Uj - 16(a fc a£ + a k + ag). 



/16 

Here {,4, B, C} = ABC + ACB + BAC + BCA + CAB + CBA. 
The eighth order functional relation is 



E 

i,j,k,l 



~^LijL k i — —{Li k , Lu, Lj k , Lji} — 7^{Lij, Li k , L k i} — —{L^, L^, L k i} 

1 /l 2 \ rr 2 r A3 3 \ r2 

+ g I 2 + ) \ L ij L ikLjkf + -LijL k i - I - - -afe - -a; - a fc aj I 









+ 








2 








+ 


a 


+ 










2 




+ 


—a, 
3 





1 \ rr r 1 I 4 4 7 

-a; I {L ik , L jk \ + I -a k + -a t + -( 
4 

afcaz + 2aiaja k + -ajC 
3 



0. 



Here {A, B,C, D} is the 24 term symmetrizer of 4 operators and the sum is taken over all 
pairwise distinct i, j, k, £. For the purposes of the representation, it is useful to redefine the 
constants as a,i = bf — \. 

We note that the algebra described above contains several copies of the algebra generated 
by the corresponding potential on the two-sphere. Namely, let us define A to be the algebra 
generated by the set {L^, 1} for all i, j = 1, . . . ,4 where / is the identity operator. Then, we 
can see that there exist subalgebras A k generated by the set {Lij, 1} for i, j ^ k and that these 
algebras are exactly those associated to the 2D analog of this system. Furthermore, if we define 

H k = ^ L ij ~ 

then H k will commute with all the elements of A k and will represent the Hamiltonian for the 
associated system. For example, take A4 to be the algebra generated by the set {L12, £13, L23, /}. 
In this algebra, we have the operator H4 = L\i + L13 + L23 + (3/4 — b\ — b\ — b^)I which is 
in the center of A4 and which is the Hamiltonian for the associated system on the two sphere 
immersed in M 3 = {(x±, X2, £3)}- 

Next we construct families of finite dimensional and infinite dimensional bounded below ir- 
reducible representations of this algebra that include those that arise from the bound states of 
the associated quantum mechanical eigenvalue problem. At the same time we will construct 
models of these representations via divided difference operators in two variables s and t. Im- 
portant tools for this construction are the results of [12] giving the representations of the A k s 
and known recurrence relations for one-variable Wilson and Racah polynomials. 




3 Review of Wilson polynomials 

Before we proceed to the model, we us present a basic overview of some of the characteristics 
of the Wilson polynomials |41] that we plan to employ in the creation of our model. The 
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polynomials are given by the expression 

) n (t 2 ) = w n (t 2 ,a,/3,7, 5) = (a + f3) n (a + 7 ) n (a + 8) r 



i". 



p ( —n, a + P + j + 5 + n — 1, a — t, a + t 
4 3 \ a + (3, a + 7, a + <5 

= (a + /3) n (a + 7 )n(a + 5)„$^'^ (t 2 ) , 

where (a) n is the Pochhammer symbol and 4^3(1) is a generalized hypergeometric function of 
unit argument. The polynomial w n (t 2 ) is symmetric in a, {3, 7, <5. 

The Wilson polynomials are eigenfunctions of a divided difference operator given as 

T *T$ n = n(n + a + p + 7 + 8 - l)* n , (3.1) 

where 

£^F(t) = F(t + A), r = ^(^-E- 1 / 2 ), 

r* = ^ [(a + + t)(7 + + *)£ 1/2 - (« - t)(/3 - 0(7 - - t)E~ 1 / 2 ] . 

See [12] for a simple derivation, 
e Wilsoi 

t 2 <$> n (t 2 ) = K(n + 1, n)$ n+ i (t 2 ) + if (n, n)$ n (i 2 ) + JC(n - 1, n)$ n _i (t 2 ) 



The Wilson polynomials <E> n (t 2 ) = $i Q '^' 7,< ^(t 2 ), satisfy the three term recurrence formula 



where 

a + /3 + j + 8 + n- 1 



if(n + l,n) 



(a + /3 + 7 + 5 + 2n - l)(a + /3 + 7 + 5 + 2n) 

x (a + /3 + n)(a + 7 + n)(a + 5 + n), (3.2) 

K( _-, v + 7 + n - !)(/? + g + n - 1)( 7 + 8 + n - 1) 

ln ' nj_ (a + /3 + 7 + 5 + 2n-2)(a + /3 + 7 + 5 + 2n-l)' ^ 

n) = a 2 - lf(n + 1, n) - K(n - 1, n). (3.4) 

This formula, together with <3?_i = 0, $0 — 1> determines the polynomials uniquely. 

We can construct other recurrence relations between Wilson polynomials of different pa- 
rameters using a family of divided difference operators L^, W ,v , fj,,v = a,/3,^,8 given in 
Appendix [A] Most importantly for the model considered below, we can construct operators 
which fix n, the degree of the polynomial and which change the parameters by integer values. 
In the model constructed below, we will want to change a and 5 by integer values and keep /3, 7 
fixed. The operators which accomplish this are given by 

L afj L aj <f>t M>S) = (« + P ~ + 7 - l^-iAT.W), 

pa/3 pQ7^(a,/3,7,5) _ (" + a + + " + 7K» + + ^ ~ + 7 + 8 - 1) 

n " (a + /3)(a + 7 ) 

We give the action on the <I>i a '^' 7 ' <5 ^(t 2 ) for simplicity. For a complete exposition on the recurrence 
relations see Appendix [A} 

Finally, the weight function of the model will be based on a two dimensional generalization 
of the weight function of the Wilson polynomials. 

For fixed a, f3, 7, 8 > (or if they occur in complex conjugate pairs with positive real 
parts) |41j . the Wilson polynomials are orthogonal with respect to the inner product 

2 

dt 

r(2tt) 



(w n ,W n >) = — I w n (~t 2 ) W n > {-t 2 ) 



r(a + it)r(p + it)r(7 + u)v(8 + it) 
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8 nn >n\{a + /3 + 7 + <5 + n-l) n 

r(a+/3+n)r(a+7+n)r(a+5+n)r(/3+7+n)r(/3+5+n)r(7+5+n) 
X r(a+/3+ 7 +<5+2n) 



. (3.5) 



When m is a nonnegative integer then a + (3 = —m < so that the above continuous 
Wilson orthogonality does not apply. The representation becomes finite dimensional and the 
orthogonality is a finite sum 

(a - 7 + l) m (a - 6 + l) m ^ (2a) fc (a + l) fc (a + /3) fc (a + 7 ) fc (a + S) k 



2a + 1 m 1 - 7 - m f— ' 1 



(2a + l) m (l - 7 - 5) m ^ (l) fc (a) fc (a - /3 + l) fc (a - 7 + l) fc (a - 5 + l) k 



x w n ((a + A;) 2 )m„/((a + /c) 2 ) = 5 nn > 



(a+/3+7+5) 



2» 



. (3.6) 



Thus, the spectrum of the multiplication operator t 2 is the set {(a + k) 2 : k = 0, . . . , m}. Now, 
we are ready to determine the model. 



4 Construction of the operators for the model 

To begin, we review some basic facts about the representation. 

The original quantum spectral problem for (|2.ip was studied in [33] from an entirely different 
point of view. It follows from this study that for the finite dimensional irreducible representations 
of the quadratic algebra the multiplicity of each energy eigenspace is (M + 2)(M + l)/2 and we 
have 



L 12 + L 13 + L 23 + L U + L 2A + L 34 = ^- \2M + J2 b j + 3J + XJ fe i | J ' 

where I is the identity operator. 

Of course, for an irreducible representation, the Hamiltonian will have to be represented by 
a constant times the identity and initially for the construction of the model, we assume 



E - 1 + h ) 1 

We will obtain the quantized values of E from the model. 

We recall that each operator Lj,- is a member of the subalgebras A k for k 7^ Thus, we 
can use the known representations of these algebras, and symmetry in the indices, to see that 
the eigenvalues of each operator will be associated with eigenfunctions 4>b,m indexed by integers 
< h < m so that 

(U^m = (-(2h + h + bj + l) 2 - i + b 2 + bfj <f> h , m , (4.2) 
(Lij + L ik + L jk )(/) Km =[\~ (2m + bi + bj +b k + 2) 2 ^j 4> Km . (4.3) 
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4.1 A basis for L 13 , L 12 + L 13 + L 23 

As described above, we seek to construct a representation of A by extending the representations 
obtained for the subalgebras Ak- The most important difference for our new representation is 
that the operator H4 = L\ 2 + L13 + L 23 + 3 / 4 — (b 2 + b 2 + b 2 ) is in the center of A4 but not A. 
Hence, it can no longer be represented as a constant. We can still use the information about its 
eigenvalues to make an informed choice for its realization. 

Restricting to bounded below irreducible representations of the quadratic algebra initially, 
we see from the representations of Ai that the possible eigenvalues of H4 are given as in (|4.3|) 
and the eigenvalues of L13 are given as in (|4.2p . 

We can begin our construction of a two-variable model for the realization of these represen- 
tations by choosing variables t and s, such that 

H A = \ -4s 2 , L 13 = -4t 2 -l + 6 2 + fe|, 

i.e., the action of these operators is multiplication by the associated transform variables. From 
the eigenvalues of the operators, we can see that the spectrum of s 2 is {(— s m ) 2 = (m + 1 + (&i + 
62 + 6 3 )/2) 2 } and the spectrum of t 2 is {tj = (I + (61 + b 2 + l)/2) 2 }. 

In this basis, the eigenfunctions di^ m for a finite dimensional representation are given by delta 
functions 



dt, m (s, t) = S(t - t £ )5(s - s m ), < £ < m < M. 
4.2 A basis for L 12 , L 12 + L 13 + L 23 

Next, we construct Lyi in the model. Let f n ,m be a basis for the model corresponding to 
simultaneous eigenvalues of L12, L\i + £13 + £23 • From the representations of Ai [U] 5 we know 
that the action of L13 on this basis is given by 

Lnf n ,m= ^2 C mU,n)f jjm , (4.4) 

j'=n,ra±l 

where 

= lJb 2 -b 2 )( h + b 2 + 2 m + 2)(b 1 + b 2 + 2b 3 + 2rn + 2) 

K ' ' 2 (2n + 6 1 + 6 2 + 2)(2n + 6i + 6 2 ) 1 3 ' K ' 

C m {n, n + l)C m (n + 1, n) = 16(n + l)(n - m)(n -b 3 - m)(n + b 2 + l)(n + b 1 + 1) 

, , , , , , n . (n + m + h + b 2 + 2)(n + m + h + b 2 + b 3 + 2) 
x (n + 6x + 6 2 + 1) (2n + - + &2 + 3)(2n + - + ^ + 2)2(2n + - + ^ - 1} . (4.6) 

We already know that the bounded below representations of Aa are intimately connected 
with the Wilson polynomials. The connection between these polynomials and the representation 
theory is the three term recurrence formula (|4,4ft for the action of L13 on an L\ 2 basis, where 
the coefficients are given by (|4.5p and ([4.6p . 

We define the operator L on the representation space of the superintegrable system by the 
action of the three term recurrence relations for the Wilson polynomials given by expansion 
coefficients (pT2l - (pT4l . i.e. 



Lf n = K(n+ l,n)f n+ x + K(n,n)f n + K(n - l,n)/ n _i. 

Note that with the choices 

61 + 63 + I a 61+63 + I 
a = to, p = , 
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7 = 2 ' 2 +^2 + ^ + 2, (4.7) 

we have a perfect match with the action of L± 3 as 

C m (n + l,n) = AK(n + l,n), C m (n-l,n) = AK(n-l,n) - - + b\ + b\. 
Thus, the action of L13 is given by 

Wn = (-4L - i + 6? + /„, 

and so we see that the action of L13 on an L12 basis is exactly the action of the variable t 2 on 
a basis of Wilson polynomials. Hence, we hypothesize that L12 takes the form of an eigenvalue 
operator for Wilson polynomials in the variable t 

L12 = -ir*r t - 2(6i + 1)(6 2 + 1) + 1/2, 

where r, r t * are given as (|3.ip with the choice of parameters as given in (|4.7p . Here the subscript i 
expresses the fact that this is a difference operator in the variable t, although the parameters 
depend on the variable s. 

The basis functions corresponding to diagonalizing H4 and Lyi can be taken, essentially, as 
the Wilson polynomials 

fn,m(t,s) = w n (t 2 ,a,/?,7,(5) 5{s - s m ), 

where s m = m + 1 + (61 + 62 + b%)/2 as above. Note that w n {t 2 ) actually depends on m (or s 2 ) 
through the parameters a, 5. Also a + 5 is independent of m. Written in terms of the variable s, 
the parameters are given by 

fr 2 + l , a h + b 3 + 1 b 1 - b 3 + 1 fc 2 + l 
« = ^— + /5 = ^ , 7 = ^ , 8 = — s. (4.8) 

Note that when s is restricted to s m , these parameters agree with (|4.7p . 

Since the t<; n are symmetric with respect to arbitrary permutations of a, ft, 7, 5, we can 
transpose a and j3 and verify that w n is a polynomial of order n in s 2 . 

4.3 A basis for £13, 2v 2 4 

For now, let us assume that we have a finite dimensional irreducible representation such that 
the simultaneous eigenspaces of L 12 , L12 + ^13 + £23 are indexed by integers n, m, respectively, 
such that < n < m < M. Each simultaneous eigenspace is one-dimensional and the total 
dimension of the representation space is (M + 1)(M + 2)/2. Now we need to determine the 
action of the operators L14, L24, L34 in the model. 

A reasonable guess of the form of the operator L24 is as a difference operator in s, since it 
commutes with L13. We hypothesize that it takes the form of an eigenvalue equation for the 
Wilson polynomials in the variable s. We require that it have eigenvalues of the form (|4.2p . 

Note that when acting on the delta basis dg m , it produces a three-term recursion relation. For 
our representation, we require that that the representation cut off at the appropriate bounds. 
That is if we write the expansion coefficients of L24 acting on d,£rn: ^ 

Ludtm = B(m,m- l)^ m _ x + B(m,m)dt m + B(m,m + l)dt m+1 , 
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we require B(m, m — l)B(m — l,m) = and B(M, M + 1)B{M + 1, M) = 0. These restrictions 
are realized in our choices of parameters, 

a = t + ^— , p = -M 1, 

7 = M + 6 4 + fcl + fe 2 2 + fc3 +2, * = -t+^±l. (4.9) 
For L24 we take 

L 24 = -4f s *f s - 2(6 2 + 1)(6 4 + 1) + i. 

Here f s is the difference operator in s where the parameters are a, /3, 7, <5. 

With the operator L24 thus defined, the unnormalized eigenfunctions of the commuting 
operators L13, L24 in the model take the form g n j. where 



£i3ff£,fc = f-(2* + 61 + £> 3 + I) 2 - \ + &? + 4) 9n,k, 
Lugtk = (~(2k + b 2 + 64 + I) 2 - 2 + fe i + b l) 5n,fc, 



where < £ < M, < k < M - £, and 

ge,k = S(t - ti)w k (s 2 ,a,(3,Zf,5), (4.10) 

with t( = £ + (61 + 63 + l)/2 as above. 

For this choice of parameters, the functions (|4.10p constitute an alternative basis for the 
representation space, consisting of polynomials in s 2 , t 2 multiplied by a delta function in s. 

4.4 Completion of the model 

In this section, we finalize the construction of our model by realizing the operator L34. The 
operator L34 must commute with L12, so we hypothesize that it is of the form 

L34 = A(s)S(L a0 L ai ) t + B{s)S~\R a PR a ~<) t + C(s)(LR) t + D(s), (4.11) 

where S u f(s,t) = f(s + u,t), A, B, C, D are rational functions of s to be determined, and the 
operators L a p, R a/3 , L, R, etc. are defined in Appendix El The subscript t denotes difference 
operators in t. (Note that r t *r = (LR)t.) The parameters are (|4.8p . Here 

L a pL ai = ^—^ (a -l + t)(a + t)((3 + t)( 7 + t)T l 

+ ±t{t- ±) {a - 1 -*)(«- *)(^ - tX-r - *) T_1 

1 -j-(a - 1 + t)(a - 1 - t){fi + 1)(7 - 1 - 1) 



+ At{t- l -) {P ~ t)[l ~ t){6 ~ 1 ~ t){S ~ t)T ^ 
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j-CS - 1 - t)( 7 + - 1 + t)(S - 1 - t) 



1 



4i(i 



- 1 + t)( 7 - t)(£ - 1 - t)(S - 1 + 1) 



2 - 



Lii = At{t+ l -) i - a + t){P + t)(l + t){6 + ^ + At{t- l -) {a ~ t)W ~ t)(7 " t){5 ~ t)T ~' 

~ 4t(t + l) (a + + ^ 7 + *)(<* + *)- 4f (/_ " MP ~ -*)(*- *)■ 

On the other hand, we can consider the action of L34 on the basis (|4.10j) . Considering L34 
primarily as an operator on s we hypothesize that it must be of the form 

L 34 = A(t)T(L^L & j)s + B{t)T-\K&Br<) a + C{t){LR) s + D(t)s 2 + E(t) + nL 12 , (4.12) 

where the difference operators are defined in Appendix [B] with subscript s denoting difference 
operators in s and K is a constant. 

Finally, we express the operator L14 as 

E + $ — 1 ) / — L12 — Ll3 — -^23 — -^24 — L34. 

U ) 

By a long and tedious computation we can verify that the 3rd order structure equations are 
satisfied if and only if E takes the values 

( 4 V 

E = - 2M + ^ h 3 '< + 3 - 1 

and the functional coefficients for L34 in (|4.11j) . (I4.12P take the following form : 
(2M + h + b 2 + 63 - 2s + 2)(2M + 61 + 6 2 + ^3 + 26 4 + 2s + 4) 



A(s) 
B(s) 



2s(2s + 1) 

(2M + 61 + 6 2 + b 3 + 2s + 2)(2M + h + 6 2 + b 3 + 26 4 - 2s + 4) 



C(s) = -2 + 



2s(2s - 1) 

2(2M + bi + 6 2 + b 3 + 3)(2M + &i + 6 2 + ^3 + 26 4 + 3) 

4s 2 - 1 



^s) = 2s^-2( 2M + &1 + b2 2 +b3 + fe4 + 4 ) - ibl + h)2 2 +bl + b Kb 3 + b i + 2M + 3 



A(t) 
B(t) 



((h + b 2 + l) 2 - 6§)(2M + 61 + b 2 + b 3 + 3)(2M + 61 + b 2 + b 3 + 2b 4 + 3) 
+ 2(4s 2 - 1) 

(61 - 63 + 2t + l)(6i + 1 + 63 + 2t) 



2t(2t + 1) 
(&i-&3-2i + l)(&i + l + & 3 -2t) 



C(t) = 2 + 



2t(2* - 1) 
2(61 " bi) 



4t 2 -l ' 

= 2, (4.13) 

and k = —4. The expression for E(t) takes the form -E(i) = fJ.i + fi 2 /(4:t 2 — 1) where Hi, \i 2 are 
constants, but we will not list it here in detail. 

For finite dimensional representations, we have the requirement that M be a positive integer 
so we obtain the quantization of the energy obtained previously (|4,1D . 
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4.5 The model and basis functions 

We shall now review what we have constructed, up to this point. We realize the algebra A by 
the following operators 




L 12 = -4t*t* - 2(6i + 1)(6 2 + 1) + ^ L 24 = -4f s *f s - 2(b 2 + 1)(6 4 + 1) + i 
L 34 = A(s)S(L a 0L a r f ) t + B(s)S^(R afS R^) t + C(s){LR) t + £>(s), 



where the parameters for the Tt operators are given in (|4.8p . the parameters for the operators f s 
are given in fj4.9|) and the functional coefficients of L34 are given in ()4.13p . The operators L23, 
L14 can be obtained through linear combinations of this basis. 

Using Maple, we have verified explicitly that this solution satisfies all of the 4th, 5th, 6th 
and 8th order structure equations. 

We have computed three sets of orthogonal basis vectors corresponding to diagonalizing three 
sets of commuting operators, {^13,^4}, {Li2,Hn} and {L13, £24}, respectively, 

dt,m(s, t) = 5(t - ti)8(s - s m ), < £ < m < M, (4.14) 
fn,m(s,t) = w n (t 2 , a, p,j,5)5(s - s m ), 0<n<m<M, (4.15) 
g e , k (s,t) = w k {s 2 ,aj,^,5)5{t-te), 0<£<k + £<M. (4.16) 

We also have a nonorthogonal basis given by 

h n>k (s, t) = t 2n s 2k , 0<n + k<M. 

Recall that the spectrum of the variables s, t is given by 

t£ = 1 H , s rn = — to + 1 H , < £ < to < M . 



We finish the construction of the model by computing normalizations for the basis f n ,m, 
and g^ m and the weight function. 

5 The weight function and normalizations 

We begin this section by determining the weight function and normalization of the basis functions 
in the finite dimensional representations. Later, we shall extend the system to the infinite 
dimensional bounded below case. 

5.1 The weight function and normalization 

of the basis di, m (s,t) = 5(t — ti)d(s — s m ) 

We consider the normalization for the di^ m = 5{t — te)5(s — s m ) basis for finite dimensional 
representations where 

, , h + 63 + 1 / , , , b 1 + b 2 + b 3 \ 

te = £ H , s rn = — [m + 1 H , < £ < to < M . 
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In order to derive these results we use the requirement that the generating operators Ljj are 
formally self- adjoint. 

Consider a weight function uj(t,s) so that 



(f(t,s),g(t,s)) = J J f(t,s)g(t,s)u)(t,s)dsdt, 

then we assume that the basis functions are orthonormal with 

{c£, m 5(t - ti)5(s - s m ),C£^ m d(t - t' e )S(s - s' m )) = 6 n 



J m,m' 



which implies that cj m uj(te, s m ) = 1. The adjoint properties of L\ 3 and L24 provide recurrence 
relations on the C£ m . That is 

{5(t -t t - l)S(s - s m ), L 13 5(t - t e )S(s - s m )) 

= (Li 3 S(t -t e - l)5{s - s m ),5(t - ti)S(s - s m )) 



implies the recurrence relation 

4+i, m = (£ + l)(l + b 3 +£)(m-e + b 2 )(m + £ + b ± + 63 + 2)(2l + b 1 + b 3 + 1) 
c l m ~ (m-£){l + b 1 + b 3 +l){m + £ + b 1 +b 2 + b 3 + 2)(2£ + bi + b 3 + 3) 

Similarly, the self-adjoint property of L 2 4 

(5(t - t t )S(s - s m + 1), L 2 4<5(* - t e )5(s - s m )) 

= (L2i5(t -t £ - l)5(s - s m + l),5(t - ti)5(s - s m )} 

implies the recurrence relation 

cj, m +i _ (M -£ + b A ){m + £ + 61 + b 3 + 2)(M + m + b 1 + b 2 + b 3 + 2) 
c lm ~ (M + m + bi + b 2 + b 3 + b 4 + 3)(m + £ + bi + b 2 + b 3 + 2) 
{m-£ + l)(2m + 2 + b 1 +b 2 + b 3 ) 
(M — m){m - £ + 1 + b 2 ){2m + 4 + b x + b 2 + b 3 ) ' 

Putting together (|5.ip and (15. 2p we obtain 



(5.1) 



(5.2) 



c jm = (1 + b 3 Ml + b A ) M {M + h + b 2 + b 3 + 3) m (2 + h + b 3 ) m+e 

c 2 0fi (1 + b 2 ) m -e( 1 + b i)e(l + h + b 3 ) e (l + b A ) M - m {M + h + b 2 + b 3 + b 4 + 3) m 
(M - m)!(m - £)!^!(2 + 61 + 6 2 + b 3 )(l + bi + ft 3 ) 
X M\{2m + 2 + b l + b 2 + b 3 ){2£ + l + b l +b 3 )(2 + b l + b 2 + b 3 ) m+ t { ' j 

which gives the value of the weight function for the spectrum of s via w{ti, 5 m ) — Cp . 
5.2 Normalization of the w n (t 2 )d(s — s m ) basis 

Next, we use the orthogonality of the Wilson polynomials to find the normalization of the fn,m 
basis in the finite dimensional representation. 

Assume the normalized basis functions have the form 

fn,m(s,t) = k n , m w n (t 2 , a, /3,7,<5) 8(s - s m ), < n < m < M. 

When evaluated at s = s m , the parameters are given by 

61 + 63 + I a 61+63 + I 
a = m, p = , 
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61 - 63 + 1 61 + 63-I , , , ^ 9 f , A , 
7 = 5 ' 2 +^2 + "i + 2, (5.4) 

and satisfy a + (3 = —m < 0. Thus, the Wilson orthogonality is realized as a finite sum over the 
weights of t 2 . However, the weight of the variable t is given by tg = I + /3 and we must adjust 
the equation for the Wilson orthogonality (|3,6p by permuting a and (3. This is allowed since 
the polynomial and the requirement a + /3 = —m are symmetric in the two parameters. In this 
form the Wilson orthogonality is given over the spectrum of the multiplication operator t 2 as 
the set {(/3 + ^) 2 : £ = 0,...,m} 

, v Cg-7 + l)mQ3-a + l)m 

)= in o 1 1 \ /1 In 5 - 5 ) 

(2p + l) m (l - 7 - 5) m 

X to ™^ - « + TO - 7 + 1M/3 - ^ + D^" ((/3 + ^ ) 
= S nn >n\{n + Q + /3 + 7 + J- l) n (q + /3) n (/3 + 7 ) n (/3 + £) n (a + 7 ) n (q + £) n ( 7 + <5) n 

(a + /3 + 7 + 5) 2n 

In light of this orthogonality, we hypothesize that the weight function is given by 



(f(t, s),g(t, s))= JJ^ f(t, s)g(t, s)w(t, s)S(t - t e )5(s - s r 
= f(te, s m )g(ti, s m )uj(ti, s m ) 



£,m 

and so we look for normalization constants so that 



(fn,m(s,t),fn',m'(s,t)) = j j ^ fn,m(s,t)f n i tm >(s,t)w(t, s)6(t - t e )5(s - S m ) 

£,m 

= 5 m ,m' k n ,Tnk n ',m w n(tj)w n >(tj)w{t£, S m )S(t - ti)S(s - S m ) 
i 

= &m,m'Sn,n'- (5-6) 

The orthogonality (|5.5p in terms of the choices of parameters (I5.4p is given by 

(2 + 6 1 + 6 2 ) 2n (m-n)!(l + 6 3 ) m—n 



~ (n + b 1 +b 2 + l) n (l + 6i)„(l + 6 2 ) n (2 + 61 + 6 2 ) m+n (2 + 61 + 6 2 + 6 3 ) m +„ 
(1 + 6i)/(l + 61 + 63)^(1 + b 2 ) m -e(2 + h + b 2 + h) m+£ 

X 



£\(m-£)\(l + 63)^(2 + 61 + 63)^ 

(2£ + 61 + 63 + I) , 2 , , 2 . /r 
x {1 + h + h) V fa 2 K^i)- (5.7) 

The weight function (|5.3p can be rewritten as 

s = M!(l + h) M - m (M + 61 + 6 2 + 63 + 64 + 3) m (2m + 2 + 61 + 6 2 + 63) 
Wl * Smj (M - m)!(l + 6 4 ) M (M + 61 + 6 2 + 63 + 3) m (2 + 61 + 6 2 + 63) 1 ' J 

< (1 + 6i)i(l + 61 + 63MI + 6 2 ) m _i(2 + 61 + 6 2 + fr 3 )m+i(2l + I + 61 + 63) 
£\(m-e)\(l + 63)^(2 + 61 + 63)^(1 + 61 + 63)40 

We can now solve the equation f)5.6[) for /c njm by comparing (|5.8p and (|5.7p to obtain 

fen,m _ (I + 64MI + 63) m—n(M + 61 + 6 2 + 6 3 + 3) m (2 + 61 + 6 2 + b%)(2 + 61 + 6 2 ) 2n 
c|~T ~ n\M\{\ + 6 4 )M-m(« + 61 + 6 2 + l) n (l + 6i)„(l + 6 2 ) n (2m + 2 + 61 + 6 2 + 63) 
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(m — n)\(M — m)\ 



(M + 6 X + 6 2 + 6 3 + 6 4 + 3) m (2 + bi + b 2 ) m+n (2 + h + b 2 + b 3 ) m+n ' 
With this normalization the basis functions f n ,m(s,t) are orthonormal. 
5.3 Normalization of the w k (s 2 )5(t — t t ) basis 

Next, we use the orthogonality of the Wilson polynomials to find the normalization of the g Utk 
basis in the finite dimensional representation. We take the normalized basis functions to be 
given by 

9£,k( s , t) = hi : kWk{s 2 , a, (3, 7, S)S(t — t/), < n < M, < k < M - £. 

Again, we want to show that there exist normalization constants ht jk so that the following 
holds: 

{k,k{s,t),g e ^{s,t)) = J J ^2g£ tk (s,t)g# )kl (s,t)w(t,s)5(t - t e )5(s - s m ) 

£,m 

= ^ hi !k hg ik m k (s 2 m )w k >(s 2 m )w{te., s m )S(t - t t )5(s - s m ) 



in 



When restricted to t = the parameters a, $, 7, 5 become 

~ B . , bi + b 2 + b 3 ~ 61 + 62 + 63 
a = £+l-\ , P = -M-l , 

71 r 1 7 1 r) 1 61+62+63 ~ 61-62 + 63 

7 = M + 6 4 + 2 + , S = -£ , (5.9) 

and so we have a + [3 = —M + £ < 0. Also, the spectrum of the variable s 2 is given by the set 
{(m+l+ fcl+j f +b3 ) 2 : m = £, . . . , M} which we can write as {((m-£)+a) 2 : m—£ = 0, . . . , M—£}. 
The Wilson orthogonality can be written in terms of the choice of parameters (|5.9p as 

s / = (M-fc + £)!(3 + 6i + 6 2 + 6 3 ) M 
(3 + 61 + 6 2 + 63 + bi) M +k-l 

(2 + 61 + 6 3 )*/+i-fc(3 + 61 + 6 2 + 63 + 6 4 ) M (2 + 6 2 + 64 + k) k 
X (1 + 6 2 ) fc (l + 6 4 ) fc (l + 6 2 + 64 + k) k (2 + 6 2 + b 4 ) M -i+k 

x ^4 (1 + 6 4 ) M -m(M + 61 + 6 2 + 63 + 64 + 3) m (2 + 61 + 6 2 + b 3 ) m+e 



(M - m)\(m - £)\{M + 61 + 6 2 + 63 + 3) m (2 + 61 + 6 3 ) r 



(1 + b 2 ) m -e(2m + 2 + 61 + 6 2 + 63) , 2 , , 2 , 
X (2 + 61 + 62 + 63) Wk M Wk ' W ' 

where the index being summed over is m = £,..., M instead of m — £ = 0, . . . , M — £. 
Comparing this orthogonality with the weight function (|5.3p written as 

, v = (1 + b 4 ) M -m(M + 61 + 6 2 + 63 + 64 + 3) m (2 + 61 + 6 2 + b 3 ) m+e 
Wl * Sm) (M - m)!(m - £)\{M + 61 + 6 2 + 6 3 ) m (2 + 61 + 6 3 ) m+ , 

(1 + 6 2 ) m _ f (2m + 2 + 61 + 62 + 63) 
(2 + 61 + 6 2 + 63) 
M!(l + 6i)^(l + 61 + b 3 )t(2£ + I + 61 + 63) 



£!(1 + 63)^(1 + 6 4 )m(2 + 61 + 6 2 + 6 3 )(1 + 61 + 6 3 )c 2 ' 



16 



E.G. Kalnins, W. Miller Jr. and S. Post 



the normalization constants are determined by the requirement 

M 

h,k> = he,kk,k ^2 S ™) w k ( S m) w k {sfn) 

m=£ 

for l^k. The proper choice of normalization is 

ht, k = (M - k + l)!(2 + bi + b 3 ) M+k+e {2 + b 2 + b 4 ) 2k 

4,o i 1 + 6 2)fc(l + &4)*(1 + b 2 + b A + k) k (3 + h + b 2 + 6 3 ) M (2 + b 2 + b A ) M +k-i 

£\(l + b 3 ) e (l + b 4 )M(l + bi + b 3 ) 

X M!(l + bi) e (l + h + b 3 ) e {M + 3 + 6i + 6 2 + 6 3 + b 4 ) m {2t + 1 + 6i + 6 3 )' 

With this normalization, the gi t k(s,t) form an orthonormal basis. 
5.4 The norm of 1 

Throughout the previous analysis, the weight and normalization have been proportional to an 
arbitrary overall factor co,o- We can fix this constant by requiring that the function 1, belonging 
to the basis h n k{s,t) of monomials in s 2 , t 2 , be normalized to length 1. We compute this by 
using the Wilson orthogonality for the Oth order Wilson polynomials (i.e. (|3.6p with k, k' = 0). 
The norm of 1 is given by 

M m 

(1, 1) = 5Z w( ^' Sm ) 
m=0 i=o 

K'f m 

= x - M!(l + b 4 ) M -m{M + h + b 2 + b 3 + b A + 3) m (2m + 2 + 6i + b 2 + 6 3 ) 
~ ha (M " m)!(1 + 54 ) M ^ M + bi+b 2 + b 3 + 3) m (2 + + b 2 + 6 3 ) 
(1 + biMl + h + ft 3 ) f (l + ft 2 )m-i(2 + fei + b 2 + fr 3 ) m+ <(2l + 1 + h + fr 3 ) 
fl(m - £)!(1 + 6 3 ),(2 + 6i + 6 3 ) m+£ (l + 6i + 63) eg, 

Evaluating the double sum gives 

(3 + 6i + 6 2 + 6 4 )m(3 + 6i + b 2 + 6 3 ) M 1 



1,1 



Thus, setting 



(1 + 6 4 )m(1 + b 3 ) M cl Q 



2 (3 + h + 6 2 + 6 4 )a/(3 + bi + 6 2 + 6 3 ) M 
C °'° (1 + 6 4 )m(1 + 6 3 ) m 

will make (1, 1) = 1. 

5.5 Infinite dimensional representations 

For infinite dimensional but bounded below representations with — m and —M + £ nonnegative 
integers and all bj real we take the inner product of two functions f(t 2 , s 2 ), g(t 2 , s 2 ) in the form 

poo roo 

(f,g)= / f(-t 2 ,-s 2 )g(-t 2 ,-s 2 )u;(t,s)dtds. 
Jo Jo 

To compute the measure ui(t, s) such that our operators Lj k are formally self-adjoint we use the 
fact that we already know the restriction of the measure to the s-constant basis (I4.15D and the 
t-constant basis ()4.16p . For consistency, we see that the weight function should be 

1 r(^+i + at + s))r(*a+i + i(-t + s ))r(^±|a±i + it )r( b -^±i + it ) 2 

(2vr) 2 r(2it) 
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r(-M - h+h+h. - i + is)T(M + h±h±h + b4 + 2 + i s ) 2 



T(2is) 



u(t, s). 



Then we can compute the norm square of the constant function f(t 2 ,s 2 ) = 1 by using ()3.5 
twice to evaluate the iterated integral: 

(1 1) = 4 r(&1 + 1)r(&2 + 1)r(&4 + 1} 



+ b 2 + b A + 3) 

x T(-M)T(-M - b 3 )T(M + h + b 2 + b 3 + b 4 + 3)T(M + h + b 2 + b 4 + 3). 

Comparing the measures that we have derived for the infinite dimensional and finite dimen- 
sional cases with the 2-variable Wilson and Racah polynomials introduced by Tratnik [38 1 139 1 HO] 
we see that they agree. Thus we have found two-variable Wilson and Racah polynomials in com- 
plete generality. 

6 Expansion coefficients 

We can easily determine the coefficients for the expansion of one of our bases in terms of another. 
Here we write the expansion coefficients in terms of the unnormalized functions. The expansion 
of the di tm basis in terms of the f n;m basis is given by 

(5(t - te)8(s - s m ),w n (t)5(s - s' m )) = 5 m , m >w n (t e )\ s=Sm Lj(ti, s m ), 
w n (t?)\ s=Sm = (-m) n (-m - b 3 ) n (l + b 2 ) n 

-n, 1 + b\ + b 2 + n, — 1 — b\ — 63 — m — £, —m + i 
-m, — m — 63, I + &2, 



X 4F3 



The expansion of the di m basis in terms of the gi^ basis is given by 

(S(t - ti)5(s - s m ),w k (s)5(t - tt)) = 5 njn u; m /_ n (s m )a;(t n ,s m ), 

w k {s m )\ t=tl = (-M + £) k (3 + h + b 2 + 63 + 64 + M + £) k (l + b 2 ) k 

—k, l + b 2 + b 4 + k, -m + i, 2 + h + b 2 + b 3 + m + 

-M + i, 3 + 6i + b 2 + b 3 + b 4 + M + £, l + b 2 , 



X 4 i ? 3 



Finally, the expansion of the f nt7n basis in terms of the gi )k basis is given by 

{w n (t)5(s - s m ),w k (s)8(t - t()) = w n (t n )\ s=Sm w k (s m )\ t= teuj(t £ , s m ). 

In order to understand the significance of these expansions in quantum theory, it is useful to 
consider the results of [43]. There the Schrodinger eigenvalue problem for the generic potential on 
the n-sphere was considered, for general n. For n = 3 it was shown that all of the eigenfunctions 
of the pairs commuting operators treated in this paper separated in some version of either 
spherical or cylindrical coordinates and were expressible as continuous multivariable orthogonal 
polynomials orthogonal on a simplex. Thus the expansion coefficients derived here represent the 
expansion of one basis of solutions of the Schrodinger eigenvalue equation in terms of another. 



6.1 A basis for L 12 , L 12 + L 14 + L 



24 



Now that we have computed the measures for our spaces of polynomials from first principles 
and established that they agree with those for the Tratnik generalization of Wilson and Racah 
polynomials to two variables |38} 139]. we can make use of known results for the Tratnik case to 
compute another ON basis for our spaces. In an appendix to [ID] the authors show that the 
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true 2-variable Racah polynomials defined by Tratnik are simultaneous eigenfunctions of two 
commuting difference operators C\, C 2 . We will identify these operators with our symmetry 
algebra and verify another eigenbasis for our representation space. 

We construct the polynomials of Tratnik [381139] and operators given in [40] via the definitions 



bi + 6 3 + 1 . h + 6 2 + 6 3 
x = 0, xx = t - , x 2 = -s-l , x 3 = M, 

Po = b 3 , Pi = 61 + 63 + I, /3 2 = 6i+& 2 + 6 3 + 2, /3 3 = 61 + 62 + 63 + 3. 

The original form of the polynomials were given in terms of the Racah polynomials which can 
be related to the Wilson polynomials via 

r n (a,b,c,d,x) = w n (a,b,c,d, (x + a) 2 ) 

with 

c+d+1 ~ c+d-l „ , d-c+1 ~ c-d+1 

a = , 6 = a , c = 6H , d= . 

2 2 2 2 

Then, the two-variable extension of the Wilson polynomials defined by Tratnik are given by 
equation (3.10) of [10] as 

R 2 (n 1 , n 2 ;Pi;x 1 ,x 2 ; M) = r m (P\ - fa — 1,02— Pi - 1, ~x 2 - 1, x 2 + Pr, x\) 

x r n2 (m + p 2 - p - 1, Ps - p 2 - 1, m - M - 1, m + & + m ; -m + x 2 ) 

with the requirements that < ri\ < n\ + n 2 < M. 

We can express the 2-variable polynomial R 2 in terms of the Wilson polynomials using the 
original parameters and variables of the model as 

R 2 (ni,n 2 ; 6;; t, s; M) = w ni (a, P, 7, 5; t 2 ) (6.1) 

/ n a + 5 a + 5 ^ r n a + S , , ^ _ a + <5 9 

x w n2 f m + /8 + — 2— ,7 + — 2~ , M + 1 + 64 + /3 + — ^— , M + 1 + /3 + — ^— , s 2 

where as in 

62 + I „ 61 + 63 + I 61-63 + I 

«= o P= n 1 7= ^ : 



In particular, note that the parameters a, 5 depend on s and so the polynomial w ni is a function 
of both s and t. 

Note that it was already demonstrated in Section 14.21 that the polynomial, w ni , is an eigen- 
function of L12. Furthermore, it is easy to see that w n2 depends only on s and so will be left 
invariant by L12 and so the 2-variable polynomial R 2 is an eigenfuction for L12. As was exhibited 
in [JOj, there is a set of two commuting difference operators whose simultaneous eigenfunctions 
are just these orthogonal polynomials. It is then natural to expect that these operators can be 
expressed in terms of the operators in our model which commute with L\ 2 , i.e., /, L12, L34, H 3 , 
and H4. 

The commuting difference operators are given as follows, via [40] . Let I{ be the operator 
which maps x% to — Xj — Pi and leaves fixed Xj for j + 1 i. Similarly, define as the operator 
which maps Xi to Xi + a and leaves fixed Xj for j ^ i. We define functions B-j ,k as 

-B- ' = Xi{Xi + Pi) + Xj + i(xj + i + Pi+i) H , 

S.f' 1 = (xj+i + Xi + - Xj + /3 i+ i - Pi), 
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Bl'° = (x i+ i - Xi)(x i+ i +Xi + ft+i), 

B]' 1 = (xn + + Xi + P i+ i)(x i+1 + Xi + p i+1 + 1), 

and further extend these functions for k = —1,0, 1 via 

B^^hiB^), B^^I i+x {B^). 

We also define by 

6° = (2xi + ft + l)(2xi + ft - 1), bj = (2 Xl + ft + l)(2xi + ft), br l = I^bj). 

Let v be some multi-index v = (v\, v<i) with v\ = —1,0, 1 and [i be a single index /i = —1,0, 1. 
Then the functions given by 

are enough to define the operators and describe the results of |40| . The operators 



r ST r w ( < a. R \ -l (A + lKft-l) \ 

^=-1,0,1 ^ ' 

£ 2 = £ a^-(xs(x s +A) + ^^i)) 

i/ i= -l,0,l ^ ' 



commute and their eigenfunctions are given by R2 (|6.ip . 
The eigenvalues of the operators C±, C2 are given by 

C1R2 = — ni(m + 61 + 62)^2, £2-^2 = -("-l + «2)(^i + n 2 + 61 + 62 + ^4)^2, 

so we hypothesize that L\ is a linear combination of L12 and the identity and £2 is a linear 
combination of #3 and the identity. In fact, it is straightforward to verify that 

L 12 = 4£ t - 26162 - 26i - 26 2 - -, 

9 

L u + £12 + ^24 = 4£ 2 - 2(6162 - 6264 - 6164) - 4(6i + 6 2 + 64) - -. 



The normalization of this basis can be found in [39] and [3D]. Thus, we have shown that Tratnik's 
version of two-variable Racah polynomials corresponds to the L12, £12 + -^14 + -^24 eigenbasis. 



7 Conclusions and discussion 

We have demonstrated explicitly the isomorphism between the quadratic algebra of the generic 
quantum superintegrable system on the 3-sphere and the quadratic algebra generated by the 
recurrence relations for two-variable Wilson and Racah polynomials, and have worked out the 
basic theory for physically relevant (boundstate) infinite as well as finite dimensional repre- 
sentations of the algebra. The 6 generators of the quadratic algebra break the degeneracy of 
the energy eigenspaces and, via various choices of commuting pairs of operators, allow one 
to describe unique bases. The eigenbases for the commuting pairs {L\2,Hi}, {£13,^4} and 
{Li2,H^} correspond to separation of the original quantum mechanical Schrodinger eigenvalue 
equation in various polyspherical coordinates, whereas the eigenbasis for {Li3,Z/24} corresponds 
to separation in cylindrical coordinates, see [33] . 
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Natural questions here are: what is the origin of these models of the symmetry algebra action 
and how can we determine when there is a differential operator model, a difference operator 
model or some other model? Clearly, the models are associated with the spectral resolutions 
of systems of commuting operators in the symmetry algebra. In [13] we showed how difference 
and differential operator models can be suggested by analysis of the corresponding classical 
systems, and these ideas are relevant here. In [44] we developed a recurrence relation approach 
for differential operators that allowed us to derive difference equation models for 2D quantum 
systems and, again, this approach should generalize to 3D quantum systems. Also, there is an 
obvious connection between the existence of models and bispectrality [40]. Another issue is that 
all models that we know of for quantum symmetry algebras of 2D and 3D superintegrable systems 
are associated with commuting operators whose simultaneous separated eigenfunctions are of 
hypergeometric type. Do there exist models with commuting operators whose simultaneous 
separated eigenfunctions are not hypergeometric? 

It is suggested by our method that most of the quadratic algebras for all Stackel equivalence 
classes of 3D second order quantum superintegrable systems on conformally flat spaces should be 
obtainable by appropriate limit processes from the quadratic algebra associated with the generic 
superintegrable system on the 3-sphere, namely that generated by the two-variable Wilson poly- 
nomials. However these limit processes are very intricate, see e.g. [33], and each equivalence class 
exhibits unique structure, so each class is important for study by itself. Moreover, within each 
class of Stackel equivalent systems the structure of the quadratic algebra remains unchanged 
but the spectral analysis of the generators for the algebra can change. We conjecture that 
this limiting process for superintegrable quantum systems is analogous to the Askey scheme for 
obtaining various families of orthogonal polynomials as limits of Askey- Wilson polynomials. 

As an example of this, in the paper [45] we studied the quadratic algebra associated with the 
quantum 3D caged isotropic oscillator. There, the Hamiltonian operator was 



and a basis for the second order constants of the motion was (with H = M\ + M<i + M3) 



We found 3 two-variable models for physically relevant irreducible representations of the quad- 
ratic algebra. One was in terms of differential operators and led to monomial eigenfunctions for 
the generators that corresponded to separation of variables in Cartesian coordinates, one was in 
terms of mixed differential-difference operators and led to one-variable dual Hahn polynomial 
eigenfunctions for the generators that corresponded to separation of variables in cylindrical 
coordinates, and the third was in terms of pure difference operators and led to one- variable 
Wilson or Racah polynomial eigenfunctions for the generators that corresponded to separation 
of variables in spherical coordinates. It can be shown that the flat space caged isotropic oscillator 
system can be obtained as a limit of the generic system on the sphere, whereas at the quadratic 
algebra level one variable dual Hahn and Wilson polynomials can be obtained as limits of two- 
variable Wilson polynomials. 

For nD nondegenerate superintegrable systems on conformally flat spaces there are 2n — 1 
functionally independent but n(n + l)/2 linearly independent generators for the quadratic al- 
gebra. It is reasonable to conjecture that the quadratic algebra of the generic potential on 
the n-sphere is uniquely associated with the (n — Invariable version of Tratnik's multivariable 
Wilson polynomials. 

Finally, these results suggest the existence of a q version of superintegrability for quantum 
systems [46] . 





M e = dj 
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A Recurrence relations for Wilson polynomials 



In addition to a three term recurrence relation, the Wilson polynomials ^^ l '^ ri ' S \t' 2 ) satisfy the 
following parameter-changing recurrence relations: 

T (a,0,j,6)^( a ,P,j,S) n(n + a + /3 + J + 5-l) ( a +i/2,/9+l/2, 7 +l/2,g+l/2) 

(a + i3)(a + y)(a + 6) n ^ 



where 



-(a,j9,7.<5) 



J_/ r l/2 _ r -l/2\ 

2y { >' 



This is a consequence of 

r[(a + y) k (a-y) k ] 
„(<x,P,ifi)ib(a,Pafi) _ 

P n 



a+ - + y 



fc-1 



fc-i. 



.( Q + /3 - 1)$ («-l/2,/3-l/2,T+W+l/2) ; 



where 



A* 



(a,/5 >7 ,«5) = J_ 

2y 



a + y-i) (/ J + y -J)^ + (a-y-i) - y - ±) ^ 



This follows from 

+ y)fc(a - 2/)fe] = + £ - fc - 1) 



1 

a y 

2 y 



Using the inner product (|3.5p . one can compute the adjoint recurrences 

r *(a+l/2,^+l/2, 7+ l/2, (5+ l/2) $ (a+l/2^ + l/2 >7+ l/2 >(5+ l/ 2 ) = (q + ^ + ^ + <j> (a ^7,*) , 

where 

r *(a+l/2,/3+l/2, 7 +l/2,<5+l/2) 

= ±- [(a + y)(/3 + y)( 7 + y)(«5 + yjT 1 ^ _ ( a _ y) (/3 _ y )( 7 - - y)^ 1 /^] 



and 



(7+l/2,<5+l/2,a-l/2,/9-l/2)^(a-l/2,/3-l/2,7+l/2,5+l/2) 

= ( ra + 7 + S)(n + a + g - 1) Q(g,i3, 7 , 5 ) 



which follows from 

( 7 +l/2,<5+l/2,a-l/2,/3-l/2) 



a ~2 +V 



1 

a y 

2 y 



= -(k + j + 6)[(a + y) k (a-y) k ] 
+ k(a + 7 + fc - l)(a + 5 + A; - l)[(a + y)jfc_i(a - 

Due to the symmetry in a, /3, 7, 5, the \x operators lead to several recurrences. Basically they 
allow one to raise any two of the parameters by 1/2 and lower the remaining two parameters 
by 1/2 in such a way that the process is essentially reversible. Moreover, composing a recurrence 
with its adjoint leads to eigenvalue equations for the Wilson polynomials. 
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B Recurrence relations for construction of the spherical 
and cylindrical models 

The spherical and cylindrical models are associated with the 8 basic raising and lowering ope- 
rators for the Wilson polynomials, as well as the three term recurrence relation. We list these 
operators here and describe their actions on the basis polynomials $ n = &i^ i '^'" f ' S \ 

R ^ = n ( n + a + P + 7 + g ~ 1) ^(a+l/2,/9+l/2,7+l/2,g+l/2) 



2. L = ^ 

2y 



(a + ft)(a + j)(a + 5) ~ n ~ x 

a - 1/2 +y)(J3- 1/2 + y)( 7 - 1/2 + y)(£ - 1/2 + y)T x / 2 
- (a - 1/2 - y)(/3 - 1/2 - y)( 7 - 1/2 - y){5 - 1/2 - y)T~ x l 2 
m n = (a + ft - l)(a + 7 - l)(a + 5 - i)^^/^-!/^-!/^-!^ 

3. L Q/3 = J- [_( a - 1/2 + y)(/3 - 1/2 + y)T x ' 2 + (a - 1/2 - y)(ft - 1/2 - y)^ 1 / 2 
2y L 



?«/9 



= — 



2y L 



5. -£^7 — — — 

2y 

L <E> = - 

6. i?« 7 = — 

2y 



.( a + /5 - 1 )$(a-l/V-l/2,TM-l/2,*+l/2)_ 

-(7 - 1/2 + y)(£ - 1/2 + ^T 1 / 2 + (7 - 1/2 - y)(* - 1/2 - y)!^ 1 / 2 

(n + 7 + <5 - 1) (re + a + /3) ^(^1/2,^+1/2,7-1/2,^-1/2) 
a + ft 

-(a - 1/2 + y)(rf- 1/2 + y)T x l 2 + (a - 1/2 - y)( 7 - 1/2 - y)T~ x / 2 

.(q, + 7 _ 1 )$(a-l/2,^+l/2,7-l/2,5+l/2)_ 

-(/3 - 1/2 + y)(<5 - 1/2 + y)T x l 2 + (ft - 1/2 - y)(<5 - 1/2 - y)T~ x l 2 
(n + ft + 5 -l)(n + a + 7) $ ( a+ i/ 2 ,M/2, 7 +i/2,w/2) 



a + 7 



a 



1/2 + - 1/2 + y)^/ 2 + (a - 1/2 - y)(<5 - 1/2 - y)T 



L a ^ n = -(a + 6- 1)^-1/^+1/2,7+1/2,5-1/2) _ 



i? Q<5 = — 



-(/3 - 1/2 + y)( 7 - 1/2 + y)T x ' 2 + (ft - 1/2 - y)( 7 - 1/2 - y)T' 

(n + ft + 7 - l)(re + a + 6) - ( a+W -i/2, 7 -i/2,a+i/a) 
a + 5 



1/2 



1/2 



Again, the three term recurrence is 

y 2 $ n (y 2 ) = K(n + 1, re)$ n+1 (y 2 ) + 2f(n, re)$ n (y 2 ) + K(n - 1, re)$ n _x (y 2 ) , (B.l) 

where the coefficients are given by (|3.2|) . (|3.3|) and (|3.4|) . 

To construct the spherical model we look for first order difference operators defining recur- 
rences for the basis functions $ n that are not just squares and that change parameters by integer 
amounts. Here are the possibilities: 

1. Eigenvalue equations. Here a, ft, 7, 5, n must be unchanged. The possibilities are 



LR$ n = n(n + a + ft + j + 5- l)$ f 
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L ap R a ^ n = { n + 1 + 8-l){n + a + ft)$ n , 
L ai R a ^ n = (n + (3 + 5-l)(n + a + 7 )$ ft , 
L a8 R aS <S> n = (n + /3 + 7 -l)(n + Q + 

2. Fix n and lower a by 1. The only possibilities are 

L^L^^'^) = (a + ft - l)(a + 7 - l^-L^+D, 
La^a**^' 7 ^ = (a + ft - l)(a + 5 - iJ^-LArH,^ 
L^L^^^ = (a + 7 - l)(a + 6 - W/r,*). 

3. Fix n and raise by a by 1. The only possibilities are 

™ ( 8 R orT $ ( a ,^ T ,J) = (n + a + ft)(ra + a + 7 )(n + ft + <5 - l)(n + 7 + £ - 1) 
n (a + ft)(a + 7 ) 

oa/J p aS at (a,P,'y,S) _ ( n + a + Z 3 ) ( n + a + + /3 + 7 - + 7 + - 1) 

^(a+l,j8,7-l,5) 

(a + P)(a + S) 

r> a jr> a S^(a,fl, 7 ,6) = ( n + a + 7) + « + g) (" + ft + 7 ~ 1) (» + ft + g ~ 1) ^(a+l,^-l, 7 ,a) 

n (a + 7 )(a + 5) 

4. Fix n and a. The only possibilities that change parameters are 

r T>a 7s> ( a ,/3,>y,5) _ (a + ft - 1) (n + a + 7) (n + ft + <5 - 1) ^ (a ^_ hl+hS) 
^ a pn v n — a + 7 n ' 

j pa**(a,A7.i) _ (a + ft - + a + £)(n + ft + 7 - 1) ^( a ,/3-i, 7 ,<5+i) 
Q/3 n a + 5 

n a + ft 

r DQ<5 a (Q,^,7,«) _ (« + 7 ~ + a + (5)(n + 7 + ft - 1) 

a + 

r TfOtpxfap^S) _ ( Q + g ~ !)( n + a + Z 3 )^ + 7 + ^-1) & (a, / 3+l, 7 ,g-l) 
a<5 n a + ft 

r Da7 ffi (a,^,7,5) _ ( a + g ~ *) ( n + a + j) ( n + Z 3 + g ~ jj ft(a,&7+l,tf-l) 

a + 7 

5. Lower n by 1 and fix a. The only possibilities that change parameters are 



n(n + q + ft + 7 + <5 - 1) ^ a ,p,-y+i,d+i) 
(a + -y)(a + 5) 
n + a + ft + 7 + J - 1) ^(afi+iaJ+i) 
(a + P){a + 5) 
n(n + q + ft + 7 + <5 - 1) ^( Q ,/3+l )7 +l,<5) 
(a + ft)(a + 7 ) 



RL a0 $ n = 5: . — . — _ L rN 



RL ai $ n = i — ^ — 



i2L Q(5 $ n = i _ . g_ L ^$ n _ x 



6. Lower n by 1 and raise a by 1. The only possibilities are 

r> 00/9(5, = n ( n + 7 + g ~ + g + ft)(n + a + ft + 7 + ^-l) (q+i^+i^/) 

(a + ft)(a + ft + l)(a + 7 )(a + 5) n ~ X 

rr«7* n(n + ft + (5 - l)(n + a + 7)(n + a + ft + 7 + (5 - 1) ( Q .+i, / 3, 7+ i,g) 

(a + 7 )(a + 7 + l)(a + ft)(a + 5) 

RR a5 $ = n ( n + ft + 7 ~ !)("- + « + ^)(^ + « + ft + 7 + ^ ~ 1) ^(a+l,/3,7/+l) 



(a + <5)(a + <5 + l)(a + ft)(a + 7) 



■n-l 
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7. Raise n by 1 and fix a. The only possibilities that change parameters are 

R a ?L<S> n = -(n + 7 + S - 2)(n + a + /3 - l)(a + 7 - l)(a + <J - l)^?' 7 " 1 '^, 
i?°TL$ n = -(n + + 5 - 2)(n + a + 7 - l)(a + /? - l)(a + 5 - 1)$&?~ 1,7, *~ 1) , 
i?^L$ n = -(n + /3 + 7 -2)(n + a + «5-l)(a + /3-l)(Q + 7- l)*^ -1 ' 7-1 ^. 

8. Raise n by 1 and lower a by 1. The only possibilities are 

L ap L<b n = -{ a + f3- l)(a + /3 - 2) (a + 7 - l)(a + 5 - l)*^ 1 '* -1,7 ^, 
L Q7 L3>„ = _( a + 7 - l)(o + 7 - 2)(a + /3 - l)(a + 5 - l)$£i 1,A7-M) , 
L a5 L<5> n = -(a + 6 - l)(a + £ - 2) (a + /3 - l)(a + 7 - l)$£r i 1A7, *~ 1) . 

In addition we can use the three term recurrence (|B. lj) and multiplication by the operator y 2 to 
both raise and lower n by 1 while fixing the other parameters. 
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